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Abstract. The shape of large on-lattice Eden clusters grown from a single seed is 
ruled by the underlying lattice anisotropy. This is reflected on the linear growth 
with time of the interface width {w ^ t), in contrast with the KPZ universality 
class {w t^^^) observed when the Eden model is grown on flat substrates. We 
propose an extended Eden model, in which the growth probability has a power law 
dependence with the number of occupied nearest neighbors. Large scale simulations 
(TV > 4 X 10^ particles) were used to determine the time evolution of w. We found that 
a suitable choice of the power exponent removes the lattice-induced cluster anisotropy 
and provides a growth exponent in very good agreement with the KPZ universality 
class. Also, the present model corroborates the results found in off-lattice simulations, 
in which the center of mass fluctuations are considered in the interface scaling analysis. 
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1. Introduction 

Since the interfaces are ubiquitous in nature, very intensive efforts are devoted to tlie 
description of tlieir pliysical properties [1]. Models wliere discrete unities are arranged in 
regular lattices constitute very useful tools due to their easy computer implementation. 
However, on-lattice growth models are frequently very sensitive to the intrinsic lattice 
anisotropy. Such a dependence is evident for clusters grown from a single particle (a 
seed) as, for example, the classical diffusion-hmited aggregation (DLA) [21 [3] and Eden 
[H [5] models. Consequently, the interface scaling analysis in these models is much less 
frequent than that for the growth from flat surfaces. However, the interface scaling in 
several experiments with radial symmetry has been considered along the past decade. 
Examples include the dynamics of several types of tumors [6l [TJ [8], the plant callus 
evolution [9], and grain-grain displacement in Hele-Shaw cells [lO], among others. In 
this scenario, the elaboration of on-lattice models without the undesirable effects of 
anisotropy becomes worthy. 

The Eden model [11] is the simplest discrete model, in which the anisotropy is 
determinant. Its original version was studied on a square lattice where the sites can 
assume two values 1 or 0, representing occupied or empty sites, respectively. The 
simulations start with a single particle at the center of the lattice and the growth 
rules are the following: at each step a peripheral site (an occupied site with at least one 
empty nearest-neighbor) is chosen at random and one of its empty nearest-neighbors 
is selected with equal probability and then occupied. Variants of the Eden rules were 
studied, and the original version is commonly called Eden B [12]. 

Despite its simplicity, the Eden model produces clusters with nontrivial interface 
scaling usually analyzed through the interface width, also called roughness. The 
roughness is defined as the root mean square deviation of the interface around its mean 
value 



where a set of distances represents the interface and f is the corresponding mean 
value. The usual behavior of the roughness is a power law growth with the time 
W ~ t^, where j3 is the growth exponent. Intensive numerical simulations of the Eden 
model grown from flat substrates provide the growth exponent /5 = 1/3 [I3l [HI [15], 
which corresponds to the Kardar-Parisi-Zhang (KPZ) universality class fT6]. However, 
the roughness of on-lattice Eden aggregates grown from a seed is ruled by the lattice 
anisotropy [11[5] and characterized by an asymptotic growth exponent (3 = 1. Zabolitzky 
and Stauffer P] simulated clusters with N ~ 10^ particles of the square lattice Eden A 
model and observed that the interface width evolves in non-trivial way. In this version, 
the growth sites (empty sites with at least one occupied nearest-neighbor) are chosen 
at random and then occupied. For small clusters, a good agreement with KPZ growth 
exponent {(3 ~ 1/3) was observed, in contrast with the linear dependence on time 




(1) 
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(/3 ^ 1) obtained for asymptotically large clusters. The value /3 = 1 is associate to the 
cluster diamond-like shape caused by the square lattice anisotropy [T7] . 

Off-lattice Eden model was simulated by Wang et al. JT8\, and the exponent 
P = 0.396, a value significantly larger than the KPZ one, was found. Recently, 
Ferreira and Alves p2j showed that the exponent reported by Wang et al. is due to the 
fluctuations of the border center of mass used as the origin for the roughness evaluation. 
Actually, the KPZ exponent is obtained when an origin fixed on the initial seed is used. 

In the present work, we analyzed large Eden clusters (A^ > 4 x 10^ particles) using 
a procedure developed in order to remove the anisotropy of DLA clusters [201 EI]- Our 
central aim is to propose a strategy for the growth of isotropic on-lattice Eden clusters 
without change their universality class, as well as to corroborate recent finds about the 
scaling of off-lattice Eden clusters. The paper outline is the following. In Section [21 the 
growth rules are presented and the results for DLA model summarized. The simulation 
results are presented and discussed in Section [3l while some conclusions are drawn in 
Section HI 

2. Neighborhood dependent growth rule 

The growth of isotropic on-lattice DLA clusters was considered by Bogoyavlenskiy |20j . 
In this rule, an empty site is chosen to be grown following the standard DLA model, 
but the growth is implemented, or not, with probability Pk given by 



where k is the number of occupied nearest-neighbors of the growth site. Although 
the preliminary simulations done by Bogoyavlenskiy suggested that the procedure was 
successful, Alves and Ferreira [21] verified that it produces patterns with diagonal 
anisotropy at the large scale simulation or the high noise reduction limits. The noise 
reduction is a tool used to enhance the cluster anisotropy and consists of a set of counters 
at each growth site that are increased by a unit every time these sites are selected for 
the growth. An empty site is occupied only after its selection for M times. Also, Alves 
and Ferreira [21] proposed a generalized algorithm with growth probabilities given by 



where n is the lattice coordination number and u an adjustable parameter. It was found 
that a suitable choice of u strongly reduces the cluster anisotropy. Indeed, the DLA 
patterns exhibit axial and diagonal anisotropies when u < and u > z/^, respectively, 
and an eight-fold structure was observed for u = z/^. In the generalized Eden model, an 
empty site is chosen to be grown using one of the Eden model versions and the growth 
is implemented, or not, following Eq. ([3]). 



(2) 




(3) 
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3. Simulations 

In this section, the square lattice case with lattice coordination n = 4 is considered. 
Figure [1] shows the borders of large Eden B clusters (linear size L = 10^) generated 
without noise reduction (M = 1) in order to illustrate the deviations from the circular 
shape. For sake of comparison, circles of diameter L centered at the seed were drawn 
(dashed line). It is evident from this figure that the border is axially or diagonally 
stretched for small and large u values, respectively, whereas an approximately symmetric 
shape is observed when Uc ~ 1. This Uc value differs from the value z/^ = 0.507 ± 0.005 
found for the noiseless limit as observed for the DLA model [21] where the crossover 
occurs at z/^ = 1.395 ± 0.005 and z/^ = 3.4 ± 0.2 for the noiseless and M = 1 cases, 
respectively. 




Figure 1. Contours of large Eden B clusters (only one quadrant is shown) without 
noise reduction for (a) = (original model), (b) i/ = 1.0, and (c) v = 2.0. In all 
patterns, the system size is i = lO'' and the total number of particles is iV ~ 10^. The 
dashed lines represent circles of diameter L used to illustrate the deviations from the 
circular shape. 

Two central points should be elucidated. The first one concerns the quantitative 
estimate for z/c, and the second one is the inquiry if the model belongs to the KPZ 
universality class for z/ = z/c. In order to analyze these points, we determined the 
interface width evolution considering the center of the lattice as the origin. Since the 
model dynamics is restricted to the border, we can do an analogy with the growth from 
flat substrates [13] and take the time proportional to the number of peripheral particles 
or, equivalently, to the mean radius f. We used square lattices containing 64000 x 64000 
sites and grew clusters with more than 4 x 10^ particles. We concentrate our analyzes 
on the Eden B model because its convergence to the asymptotic scaling is faster than 
that of the Eden A model [I2]. For the last one, ~ 10^ particles were necessary to 
enhance the diamond-like scaling [w ~ t) [4J. In Figure [2](a), the interface width is 
plotted as a function of the mean radius for z/ varying around z/ = 1. As one can see, 
the double-logarithm curves exhibit upward curvatures for both u < 1.0 and z/ > 1.0 
(inset) while a power law is observed for u ^ 1. Also, the slopes are very close to the 
KPZ growth exponent (3=1/3 when u ^ 1. 
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Figure 2. (a) Interface width as a function of the cluster mean radius for distinct v 
values and (b) mean curvature as a function oi A = ly — I'c- In (a), the main figure 
shows the curves for < I while the inset shows those for > 1. The dashed lines 
represent power law w ^ . In (b), the solid lines are the power laws R ~ (— A)°-^^ 
and R ~ /\0.78 f^j. /\ < g and A > 0, respectively. In all data, the averages were done 
over 120 independent samples. 



In order to find the null curvature point, which represents the absence of anisotropy, 
the double- logarithm data are fitted by cubic polynomials Psi^r) and the mean curvature 
is defined by 

I /"logio 

R = / K{r)dr, (4) 



where rj and r/ are the lower and upper cutoffs of the data used in the cubic fits, 
respectively, and /t(r) is the local curvature defined by 

P" 

The transient behavior for very small clusters was avoided using fit ranges from ^ 46 
to Tf ^ 3 X 10"^. The mean curvature as a function of u is plotted in figure [2]^b). This 
figure suggests a abrupt variation around the minimum of K. The minimum of h occurs 
between u = 1.001 and u = 0.999 providing the estimates = 1.000 ± 0.001 and 
f3 = 0.325 ± 0.017. Here, the fits were done for f > 100 and the uncertainty represents 
the local slope fluctuations along the distinct parts of the curves. This growth exponent 
is in very good agreement with the KPZ universality class expected for the original Eden 
model. Finally, these data strongly indicate that the curvature vanishes as a power law 
as 1/ — > Uc- 

A central feature related to the interface scaling of round clusters is the choice of 
the origin, as recently demonstrated for the off-lattice Eden model [19]. In particular. 
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Figure 3. Interface width (two methods) and CM mean distance from the initial seed 
(inset) as functions of the mean radius for v — I. The straight Unes are power law fits. 
The averages were done over 120 independent samples. 

if the width is evaluated in relation to the border center of mass (CM), the growth 
exponent is jScM — 2/5 [19], which differs from the KPZ exponent observed when the 
initial seed is used (a fixed origin). This difference was associated to the CM fluctuations 
growing faster than the interface ones. We applied this analysis to the present model 
with V = Vc 'ia order to assess the agreement with the scaling exponents obtained for the 
off-lattice model and to extend the off-lattice simulations limited to 3 x 10^ particles to 
more than 4 x 10^ particles. The results are show in Figure [3], where the interface width 
evaluated in relation to the CM {Wcm) is compared with that in relation to the seed 
{Wq). It is evident that the corresponding growth exponents are distinct. Indeed, the 
exponents are /3o = 0.324±0.017 and (3cm = 0.392±0.016. These results are in excellent 
agreement with the off-lattice Eden model [21] , reinforcing the equivalence between the 
universality classes of the models. The CM mean distance from the seed exhibits two 
power law regimes, namely Rcm ~ ^ "'^^ and Rcm ~ ^ for small and large system 
sizes, respectively. Again, the asymptotic scaling agrees with the off-lattice simulations. 

The most used version of the Eden model is that known as Eden A [12]. In this 
version a growth site (an empty site neighboring an occupied one) is chosen at random 
and then occupied. In the noiseless limit, the transition from axial to diagonal anisotropy 
is the same observed for the Eden B model, including the interface width minimum 
occurring at z/^ ~ 0.51. Without noise reduction, the criteria of null curvature gives 
Uc = 1.72 ±0.01 and growth exponents in very good agreement with the off-lattice Eden 
model. The lower precision of the parameter in comparison with the Eden B version, 
is due to the slower convergence to the asymptotic exponents and the larger Uc value. 
The last implies in a higher frequency of rejected tentatives along the simulations. 
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Figure 4. Differences between diagonal and axial growth. The dashed vertical lines 
indicates the I'c values for Eden A (right) and B (left) models. The averages were 
done over 1000 samples of size L — 400. It is important to stress that the histograms 
remains essentially unchanged as the system size is varied. 

The difference between the z/^ values for Eden A and B models can be understood 
analyzing the growth site neighborhood at z/ ?a z/^. Lets Qk be the fraction of growth 
sites with k occupied neighbors. Considering the square lattice anisotropy, one can easily 
conclude that the growth sites with k = 1 and k = 3 produce axial growth, those with 
k = 2 produce diagonal growth, while those with k = 4 produce an inward growth which 
does not contribute for the overall anisotropy. One can suppose that if the axial growth 
probability exceeds the axial one, the anisotropy rotation emerges. So, we analyzed the 
difference between them using the relation 

Figure H] shows the curves by against v. As one can see, the v = values for which 
the diagonal and axial growth probabilities are equal [by = 0) are slightly larger than 
i>c. Indeed, we found z/* — z/g ~ 0.02 and z/* — z/g ~ 0.04 for the Eden A and B, 
respectively. Notice that these differences are small, but neatly distinguishable since the 
respective deviations from null curvature are evident as shown in Figure [2]J^b). Thus, 
we conclude that the isotropic growth of Eden clusters cannot be explained only by the 
counterbalance between diagonal and axial growths. 

4. Conclusions 

In this work, we developed a strategy to grow on-lattice Eden clusters without the 
undesirable anisotropy effects. Asymptotically large clusters generated with the original 
Eden model are ruled by a diamond-like and, consequently, a growth exponent /3 = 1, in 
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disagreement with the KPZ universahty class {(3 = 1/3) observed for the growth from a 
flat substrate. We proposed an extended Eden model, in which the growth probability 
is proportional to a power of the number of occupied nearest neighbors, Pk oc k'^, in 
analogy with a procedure developed to remove the anisotropy of DLA clusters [201 l2T] . 
We found that the choices z/ = 1.72 ± 0.01 and z/ = 1.000 ± 0.001 lead to isotropic 
patterns for the Eden A and B models, respectively, preserving the KPZ universality 
class. Also, the extended model agrees with off-lattice simulations, in which the center 
of mass fluctuations are considered [T9] . 

The present strategy constitutes a useful tool for modeling experiments with radial 
symmetry since lattice models can be easily generalized. A very important example is 
the recent report concerning the dynamics of tumors |6l[7]. These experiments indicate 
a universal dynamics characterized by the Mullins-Herring universality class [22] for 
distinct types of tumors. However, this idea is subject of recent controversies [23] . 
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